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We investigate the critical behavior of the gauge glass model for the vortex glass transition in three- 
dimensional superconductors, including screening of the interaction between vortices. A Monte Carlo 
study of the linear resistivity and a scaling analysis of current-voltage characteristics indicates that 
screening destroys the finite-temperature transition found earlier when screening was neglected. The 
correlation length exponent at the resulting zero temperature transition is found to be v = 1.05±0.1. 
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The question of whether type-II superconductors in 
the mixed state exhibit a finite linear resistance, py ln , in 
the limit of a vanishingly small current has been under 
extensive consideration in recent years. While an applied 
current perpendicular to the field generates a Lorentz 
force and hence leads to motion of flux lines producing a 
voltagea, it has been proposed by Fishem that disorder 
(point defects) collectively pins the vortices, resulting in 
a phase with vanishing linear resistance. Since disorder 
also destroys the -triangular vortex lattice predicted by 
mean field theory!., this phase is called the vortex glass. 

Much of the theoretical work has been devoted to the 
study of a simplified model, called the gauge glass, which 
has the necessary ingredients of randomness and frustra- 
tion and is believed to be in the same universality claps 
as the vortex glass. In two dimensions, experimentsu, 
and simulations.]" agree that the vortex glass transition 
only occurs at T c = 0. In d = 3 the situation isj-some- 
what less clear. Several experiments on YBCOtj have 
found convincing evidence for a finite T c phase transition, 
presumably to a vortex glass. Furthermore, simulations 
support a transition at T c > 0, although the possibility 
that T c = and the lower critical dimension is precisely 
d = 3 cannot be fully ruled outEfQ. 

Most of the theoretical work on the vortex glass has 
neglected screening of the interaction between vortices, 
which arises from coupling to the fluctuating magnetic 
field. However, sufficiently close to T c the screening 
length, A, and the correlation length, £, become com- 
parable and the system-crosses over into a region where 
screening is importantLni3. A recent domain wall renor- 
malization study of the d = 3 gauge glass model by one of 
ua. (henceforth referred to as BY) concluded that screen- 
ing apparently destroys the finite temperature transi- 
tion. For experimental systems this effect would lead to 
a rounding out of the transition close to the nominal T c , 
and the linear resistance would not vanish completely but 
become very small. So far, experiments have not had the 
sensitivity to see this effect. However, the calculations 
of BY were on very small lattice sizes, L < 4, and so 
their conclusions can only be tentative. In this paper, we 
investigate the role of screening by Monte Carlo simula- 
tions, for which much larger sizes are possible, L < 12. 
Our results confirm the conclusion of BY that screening 
destroys the finite-T transition. We are also able to ob- 
tain, for the first time, values for the exponents at the 
resulting T — transition. 



The model studied here is the gauge glass with a fluc- 
tuating vector potential 

H = -j)T] cos((/) 4 - <f>j - - Aq 1 a ij ) 
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where J is the interaction strength, (pi is the phase of 
the condensate on site i, and the sum is over all near- 
est neighbors on a simple cubic lattice. The ef- 
fects of the external field and disorder are represented by 
quenched vector potentials , which we take to be uni- 
formly distributed on the interval [0, 2ir]. The fluctuating 
vector potentials are integrated over from — oo to oo, 
with a gauge constraint V • a = 0, and Ao is the bare 
screening length. In the last term, which describes the 
magnetic energy, the sum is over all elementary squares 
on the lattice and the curl is the directed sum of the 
vector potentials round the square. 

We consider the strong screening limit, A — > 0, which 
is technically easier to simulate in the vortex represen- 
tation. This is done by replacing the cosine in Eq. (|l|) 
with the periodic Villaip_.fpDjC.tion and performing fairly 
standard manipulations___T__l_] to obtain 
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Here the are vortex variables which sit on the links 
of the dual lattice (which is also a simple cubic lattice) , 
G(i — j) is the screened lattice Green's function 



G(i-j) = J 



L 3 



1 - cxp[i k- (i-j - rA] 
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(with d = 3), and the hi are given by (l/27r) times the di- 
rected sum of the quenched vector potential on the orig- 
inal lattice surrounding the link on the dual lattice on 
which hi lies. Due to periodic boundary conditions, we 
have the global constraints Yli — _Ci n * = 0- There 
are also the local constraints, [V-n]j = [V-b]j = 0, where 
the latter just follows trivially from the definition of 
as a lattice curl. 

Clearly G(0) = and, in the limit A -> 0, G{r) = 
J(2ir\o) 2 for r ^ with corrections which are expo- 
nentially small, i.e. of order exp(— r/Ao). Because 
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^(bj — rij) = we can always add a constant to G(r) 
for all r without affecting the results. We therefore add 
— (27rAo) 2 J, as a result of which the only interaction is 
on-site, and then work in units where (2-7rAo) 2 J is unity. 
The resulting Hamiltonian then has the very simple form 
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which is convenient for simulations. Note, however, that 
TCv is not trivial because the local constraint [V ■ n] j =0 
effectively generates interactions between the n.;. Note 
also that without disorder (i.e. all b^ = 0) this-model 
is just a dual representation of the XF-modefcjtJ (with 
the Villain potential) , in which the temperature scale has 
been inverted. 

We simulate the Hamiltonian in Eq. (^) on simple cu- 
bic lattices with N = L 3 sites where 4 < L < 12. Pe- 
riodic boundary conditions are imposed. We start with 
configurations with all = 0, which clearly satisfies the 
constraints, and a Monte Carlo move consists of trying 
to create a loop of four vortices around a square. This 
trial state is accepted with probability I/(I+exp(/3Ai?)), 
where AE is the change of energy and j3 = 1/T. Each 
time a loop is formed it generates a voltage AQ = ±f 
perpendicular to it's plane, the sign depending oa the 
orientation of the loop. This leads to a net voltage 



V(t) = —I v (t) with I v (t) 
2e 
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to simulate a non-uniform system, so instead we effec- 
tively assume that the current is the same everywhere so 
each vortex feels a Lorentz force x J. The scaling be- 
havior of the response to such a perturbation should be 
the same as that derived earlier for response to an actual 
transport currentE£l. We can therefore use this approach 
to determine critical exponents, which is our objective. 
The Lorentz force biases the moves and sets up a net 
flow of vortices perpendicular to the current, whose time 
average gives the voltage according ta3 Eq. (||). 

To analyze our data we need to understand the scal- 
ing behavior of the T-T^-curves near a second order phase 
transition. We quickly review the scaling thcorytj assum- 
ing T c = so the correlation length diverges as £ ~ T~ v . 
The vector potential enters the Hamiltonian in Eq. (|l|) 
in the dimensionless form j4y = J A(r) • dr so A scales 
as The electric field is given by E = —dtA and so 
scales as l/(£r) where r is the relaxation time. J • E is 
the energy dissipated per unit volume and unit time and 
scales adia T/(£<V). Hence J scales like T/£ d_1 . Com- 
bining these results leads to 
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where g is a scaling function. In three dimension this 
becomes 



(9) 



where I is the vortex-current and t denotes Monte 
Carlo "time" incremented by At for each attempted 
Monte Carlo move. We will work in units where h/(2e) = 
1, and we set At = 1/3N so that an attempt is made to 
create or destroy one vortex loop per square in each di- 
rection, on average, per unit time. 

The linear resistivity can be calculated from the volt- 
age fluctuations via the Kubo formulalla 
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Near a second order phase transition the linear resistivity 
obeys the scaling lawliO 



1 (T,L)=L-( 2 - d +^p(L 1 ^(T-T c )), 



(7) 



where p is a scaling function and z is the dynamical expo- 
nent. At the critical temperature, p becomes a constant 
and therefore pn n (T c ,L) - L -(2-d+z)^ Jf we plot the 
ratio ln[/9ii n (L)/ j oii n (L')]/ln[L/i / ] against T, then at the 
point (T c , d — 2 — z) all curves for different pairs (L, L') 
should intersect and one can read off the values of T c and 
z. We will refer to this kind of data plot as the "inter- 
section method" . 

In addition to pi; n , we also measure the voltage gener- 
ated by a finite external current. In real superconductors, 
transport currents generate a non-uniform magnetic field 
because of Ampere's law, V x B = J. It is inconvenient 



From this equation we can see that the current scale, 
Jnl, at which nonlinear behavior sets in varies with T 
as Jnl ~ T 1+2v . Since the linear resistivity is defined by 



and g(Q) can be taken to be unity, we can write 
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Furthermore, we expect that near the T = transition, 
long time dynamics will be governed by activation over 
barriers. Hence we expect 

T 1+ > lin = - = Acxp(-AE(T)/T), (12) 

T 

where AE is the typical barrier that a vortex has to cross 
to move a distance £. One can define a barrier height 
exponent ip by AE ~ ~ T~^ v in terms of which 
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We are able to obtain a rough estimate for ip from our 
data of the linear resistivity. 

In a finite system, the I-V characteristics will also de- 
pend on the ratio L/£. One can generalize the scaling 
function, Eq. (|TT|), to account for finite size effects as 
follows: 
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FIG. 1. Plot of \n[ Plin (L)/p hn (L')]/\n[L/L'] versus T for 
the pure system, whose Hamiltonian is given by Eq. (^]) with 
hi — 0. This is a dual representation of the pure XY model 
(with the Villain potential) with inversion of the temperature 
scale. The curves intersect at T c — 0.331 ±0.002. At the inter- 
section point, the y-vahie is approximately —2 corresponding 
to 2 ~ 3 according to Eq. (T7T) - 
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Now we are left with a rather complicated scaling func- 
tion since it depends on two variables. We estimated v by 
determining the current where E/(Jpu n ) = 2, at which 
point non-linear effects start to become significant. De- 
noting these values of J by Jnl, then, from Eq. (|j"4|) , it 
follows that Jnl/T 1 ^ is a function of L X I V T. Collaps- 
ing the data in the appropriate plot then gives an esti- 
mate of v. We then collected data for sizes and tempera- 
tures such that L X I V T is constant so the scaling function 
in Eq. (|l4|) then depends only on one variable. 

Let us now move to the analysis of our data. In Fig. |l| 
we show data of the linear resistivity for the pure case 
plotted according to the intersection method. All curves 
intersect at about T — 0.331 ± 0.002 which is pin— ex- 
cellent agreement with the well established valuetHrliil of 
T c = 0.33. The y-axis value at the intersection point 
is 1 — z according to Eq. (|t]), from which we estimate 
z ~ 3. This is larger than one would expect for a model 
with relaxation dynamics where usually z ~ 2. Perhaps 
corrections to finite size scaling make the estimates of T c 
and z slightly inaccurate, with the effect on z being more 
pronounced since the slopes of the curves in Fig. ^ are 
quite steep. Nonetheless, we are interested in whether 
or not a transition occurs rather than precise values of 
exponents, and the method seems to be a reliable in pre- 
dicting this, at least for the pure case. 

For comparison, Fig. |^ shows the intersection method 



FIG. 2. Plot of ln[pi in (L)/pi in (L')]/ln[i/i'] against T for 
the gauge glass model whose Hamiltonian is given by Eq. (Q) . 
In contrast to the data for the pure system in Fig. hi there is 
no intersection over the entire temperature range, indicating 
the absence of a phase transition. 



applied to the gauge glass model. As one can see, there 
is no apparent intersection over the entire temperature 
range that we have been able to simulate, i.e. down to 
T = 0.1 for L = 12 and down to T = 0.07 for L < 8. This 
rules out a transition down to about 1/5 of the critical 
temperature of the pure system, and strongly suggests 
that the phase transition is at zero temperature. 

Fig. H shows a scaling plot of the nonlinear I-V char- 
acteristics according to Eq. (|l4]). From the scaling of Jnl 
we estimated v to be approximately 1.1, and then chose 
temperatures and sizes for the data in Fig. || such that 
L X I V T, and hence the second argument in Eq. (|l4|), re- 
mained roughly constant. The data is seen to scale very 
well. Combining this data together with scaling plots for 
other sizes and temperatures we obtained the overall es- 
timate v = 1.05 ± 0.1. We also tried to scale our data 
with an appropriate scaling function for finite T c and 
found that scaling works satisfactorily with T c ~ 0.04 
and v ~ 1.1. Therefore, we conclude that the transition 
is very likely to occur at T c = 0, although a finite but 
extremely small T c cannot be completely ruled out. 

Finally, Fig. || shows data for T 2 pu n versus 1/T in 
a log-linear plot according to Eq. (|l3|). The data for 
L = 12 follows almost a straight line, i.e. shows close 
to Arrhenius behavior, which corresponds to a barrier 
exponent -0 = according to Eq. jl^). It is therefore 
possible that the barriers only increase logarithmically 
as the temperature approaches zero. However, such be- 
havior is difficult to observe in finite size systems over a 
modest range of temperatures. It is also possible that we 
are only measuring an effective exponent and that the 
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FIG. 3. Scaling plot of the /-^-characteristics with T c = 
and v = 1.05, according to Eq. jl4|, choosing sizes and tem- 
peratures such that Z/ 1/ 'T is roughly constant. 

true ip is larger than zero. 

In summary, we have performed a Monte Carlo study 
of the d = 3 gauge glass model in the vortex representa- 
tion with strong screening. We have analyzed the linear 
resistivity and I-V characteristics by means of finite size 
scaling. There appears to be no transition of this model 
into a vortex glass state at finite temperature — consis- 
tent with the findings of BY — and the correlation length 
exponent at the resulting zero temperature transition is 
found to be v — 1.05 ± 0.1. However, the possibility 
that T c is small but non-zero cannot be completely ruled 
out since we are not able to equilibrate a range of system 
sizes below about 1/5 of the transition temperature of the 
pure system. Assuming T c — 0, the presumed finite tem- 
perature vortex glass transition in the three-dimensional 
gauge glass is destroyed by screening effects. This means 
that the linear resistance in experimental systems would 
not strictly vanish, though it would become extremely 
small for T near the apparent T c . It would be interesting 
to look for this effect experimentally. 
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